Critical behavior of the 5 = 3/2 antiferromagnetic Heisenberg chain 
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Using the density-matrix renormalization-group technique we study the long- wavelength properties 
of the spin 5 = 3/2 nearest-neighbor Heisenberg chain. We obtain an accurate value for the spin 
velocity v = 3.87 ± 0.02, in agreement with experiment. Our results show conclusively that the 
model belongs to the same universality class as the 5 = 1/2 Heisenberg chain, with a conformal 
central charge c = 1 and critical exponent rj = 1. 



The study of one-dimensional spin-S Heisenberg 
Hamiltonians has received considerable attention in re- 
cent years. Their properties are important to experimen- 
talists because they describe a number of materials with 
magnetic ions arranged in chains ffijM. The first quasi-lD 
antiferromagnets with spin 3/2 which found experimental 
interest are CsVCl 3 f|] and AgCrP 2 S 6 |J]. For theorists 
the spin models are among the simplest Hamiltonians 
where quantum fluctuations are crucial. However, only 
the special case S=l/2 can be solved exactly ||. In this 
case Luther and Peschel || found that the spectrum is 
gapless, there is no magnetic ordering in the ground state, 
and the spin correlation function ui(l) decays according 
to a power law \u)(l)\ ~ l/l v with r) = 1. Later Haldane 
proposed that for half-odd-integer (h.o.i.) spin the 
Heisenberg Hamiltonian would behave as in the S=l/2 
case, i.e., no energy gap and power-law decay of the spin 
correlation function, while in the integer case the spec- 
trum would have an energy gap and the spin correla- 
tions would decay exponentially. Early numerical work 
confirmed this conjecture and power-law decay of 

the spin correlation was found for half-odd-integer spin, 
yet the exponent rj seemed to be S dependent. When 
the powerful tools of conformal invariance started to be 
applied to one-dimensional quantum systems it became 
clear that the integrable S=l/2 Hamiltonian at low ener- 
gies is equivalent to the SU(2) symmetric Wess-Zumino- 
Witten (WZW) model with topological coupling constant 
k = 1 [pLCfl . In this model, k is related to the central 
charge c by c = 3k/ (2 + k) and the critical exponent rj 
is given by 3/(2 + k). It was also found |l(J that a new 
class of integrable antiferromagnetic Hamiltonians with 
arbitrary spin |llj were equivalent to the WZW model 
with k =2S. 

Based on the SU(2) symmetry of the non-integrable 
h.o.i. spin Hamiltonians Schulz suggested the fol- 
lowing possible scenarios: 

(a) The h.o.i. spin-S Heisenberg Hamiltonian belongs to 
the same universality class as the spin-S integrable model. 
In this case both models would have c = 3k/ (2 + k) and 



rj = 3/(2 + k) and k = 25. This scenario was favored 
initially by Affleck Q. 

(b) The spin-S integrable models represent unstable mul- 
ticritical points in the phase diagram and the stable criti- 
cal point is defined by the S=l/2 Hamiltonian. Thus, the 
h.o.i. spin-S Heisenberg model is expected to be equiva- 
lent to the k = 1 WZW model having c = 1 and rj = 1, 
independent of the value of S. This position was favored 
by Schulz and later on by Haldane and Affleck |l^ , [l5| . 

Since the derivations of scenarios a) and b) involve 
several assumptions and approximations, unbiased nu- 
merical studies are needed to decide which one is correct. 
The S=3/2 Heisenberg Hamiltonian is the simplest of the 
higher half-odd-integer spin models. Since there are 4 de- 
grees of freedom per site the longest chains that can be 
exactly diagonalized with present computer capabilities 
have 14 sites. Previous numerical studies conducted in 
such small chains gave conflicting results [ p~6|JTr| . Moreo 
p6[ , analyzing the behavior of spin correlation functions 
in lattices with N < 12 found rj s» 0.6, while Ziman 
and Schulz JTtJ obtained rj w 1 by studying the behavior 
of the lowest triplet and singlet state energy gaps also 
in chains with N < 12. A Monte Carlo calculation of 
the spin correlation function |l8j in chains with up to 
128 sites provided exponents in disagreement with the 
SU(2) WZW scaling theory. The analysis of small sys- 
tems is hampered by the presence of logarithmic correc- 
tions which modify the power-law dependence of various 
quantities in a significant way. In this paper we will elu- 
cidate this problem numerically by studying sufficiently 
large systems that cannot be solved with exact diagonal- 
ization, and by Monte Carlo techniques only with con- 
siderable statistical error. 

For this purpose we will use the density-matrix 
renormalization-group technique (DMRG) |l9| keeping 
up to rn = 1200 states per block and considering chains 
up to N = 60 sites. For the largest systems considered 
we estimate that the ground state energy has an absolute 
error of less than 7 x 10 ~ 6 , by considering different m val- 
ues. Our results for small systems, N < 14, agree with 
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exact calculations [^0|,|T],[| and the ground state energy 
is exact up to N = 12. 

To calculate the conformal central charge c (or equiva- 
lently k) we account for the finite-size corrections of the 
ground-state energy per site, Eq/N, as derived from the 
WZW-theory g§ 
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The coefficient B depends on the effective coupling con- 
stant (see below) and can be dropped only in the large-N 
limit [£2| . The last higher-order term reflects a finite size 
correction to the conformal charge. In principle Eq.(l) 
may be fitted with 5 free parameters, however we found 
that more accurate results are obtained by an indepen- 
dent calculation of the spin velocity v and of B. 

The spin velocity can be determined numerically 
from the excitation energy A q = E(2n/N) — E(0) = 
v sin(27r/iV), where E(2tt/N) is the energy of the first ex- 
cited state with the lowest non-zero momentum obtained 
from the spin excitation spectrum p3| and the Davidson 
algorithm |Q (with estimated absolute error ~ 10 -2 ). 
Studying the behavior of A 9 /iV vs sin(27r/iV)/7V, as 
shown in Fig. 1, to minimize the errors, we obtain a 
rather precise value for the spin velocity v = 3.87 ± 0.02. 
This implies an enhancement factor of 1.29 with respect 
to the spin-wave result v sw = 2S = 3. This factor is sig- 
nificantly lower than the value ir/2 for the spin- 1/2 chain 
(where v sw = 1) thus implying that the quantum renor- 
malization of the low-lying spin excitations is weaker. 
The velocity is in good agreement with the measured 
spin-wave velocity v — l.26v sw for the spin-3/2 chain 
system CsVCls ||. The value reported for AgCrP 2 SQ, 
v = 1.53v su , Q , appears too high in view of our results. 
Our v is higher than the value v = 3.64 ± 0.08 given 
by Yamamoto [g5| . The latter was calculated using a 
world-line Monte Carlo approach with larger error bars 
in the excitation energies (~ 10% as estimated from the 
figures). 

The value of the coupling constant g can be determined 
directly by writing Eq. 1 in its equivalent form: 
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From Eqs. (2) and (3) and the velocity given above we 
find eoo = -2.82833(1), k = 1.00(1) and the coupling 
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FIG. 1. Excitation energy A q /N vs. sin(27r/iV)/A^. 
The linear interpolation yields the spin velocity (slope) 
v = 3.87 ±0.02 (full line). 



constant go = g(lo = 20) = 0.055(5). The error was esti- 
mated using different subsets of the data. The constant 
S depends on the smallest system included in the fit, e.g. 
for TV > 12 we found S = 20. From Eq. |, B = 0.6 and 
9eff = 0.4. 

Using the previous values of v and B and from a three 
parameter fit of Eq. |l| using N > 8 we obtain consistently 
= -2.82833, k = 1 ± 0.01 and S = 15 (Fig. 2). 
Our ground-state energy is considerably lower than the 
energy eoo = —2.8248 obtained in a recent Monte Carlo 
study Ref. p6| , where a l/A r2 -scaling of the energy was 
assumed. From our results for k we conclude that c = 
1.00 ±0.01 

To obtain the critical exponent we calculated the spin 
correlation function keeping 1000 states per block. The 
absolute error for the largest system is lower than 10~ 3 . 
We found that the scaling behavior obeyed by the spin 
1/2 chain Q| is absent in this case and we extrapolated 
the finite-size data with a quadratic polynomial in l/N 2 . 
The correlation function in the thermodynamic limit is 
shown in Fig. 3. We interpolated this correlation function 
using the theoretical prediction (l[| 



ui(l, oo ) = a 



(5) 



and obtained a = 2.06, B = 0.56 and 77 = 1.02. This 
fit is presented in Fig. 3 where we added the correlations 
for S = 1/2 to show that in both cases they have the 
same multiplicative logarithmic corrections, the coupling 
constant being smaller for S — 1/2 p2j. We note that 
a = to(la, 00) J g e f JIq and B follows from Eq. ^. Using the 
value for g e f f given above and the extrapolated value for 
u>(Iq,oo) we obtain a = 2.07 and B = 0.6. Thus the 
parameters obtained from the correlation function are 
in complete agreement, whithin error bars, with those 
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FIG. 2. Scaling of the ground state energy Eq/N. In- 
terpolation including logarithmic corrections yields the cen- 
tral charge c = 1.00 ± 0.01 and an extrapolated value 
= -2.82833 (full line). 



obtained from the scaling of the ground state energy. 

On the other hand, if we assume 77 = 1, we can analyze 
lu>(l, oo)/Zoo->(Zo, oo) = y/g{lo)/g(l)- This ratio directly 
reflects the dependence on the coupling function (22) and 
a single parameter fit yields g(lo = 20) = 0.052. 

The critical exponent 77 can also be calculated from the 
scaling of singlet and triplet excitation energies pa] : 
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Here rji/2 (i = s,t) is the scaling dimension of the cor- 
responding primary field, and Sl and Sr are the spin 
operators related to the energy levels of the conformal 
field theory describing the fixed point, with spin quantum 
numbers sl and sr. According to the WZW mapping 
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s{s + 1) - s L {s L + 1) - s r {sr + 1)], (7) 



where s is the total spin number. The triplet state s = 1 
and the first singlet excitation s = (for both cases 
sl = sr = 1/2) become degenerate in the thermody- 
namic limit. Therefore, if we take the appropriate linear 
combination of Eq. [6] for these excitations the leading 
logarithmic corrections are cancelled out: 



A = ±[(E S -E ) 
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This quantity has been used before in Ref. |17j to esti- 
mate 77 using three system sizes, i.e. N = 8, 10 and 12 
sites with the result 77 ~ 1. We reached an excellent in- 
terpolation of our data (including a S' /N 4 term), shown 
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FIG. 3. Extrapolated and normalized spin correlation func- 
tions. The fit which includes the multiplicative logarithmic 
correction gives 77 = 1.02 (full line) and coefficients that are 
in agreement with those obtained from the scaling of the en- 
ergy. The correlation function and fit of the spin— 1/2 chain 
p2| is shown for comparison. Inset: Normalized correlation 
functions multiplied by I, showing the logarithmic correction. 



in Fig. 4, leading to 77 = 1.018, b = b s - 3b t — —5.4 and 
S' = 80. 

We also obtained a good fit for the triplet finite-size 
gap (A t ), also shown in Fig. 4, but we were forced to 
add the subleading terms (~ (In \n(BN)) 2 / hi (BN) and 
S" /N A ) to obtain convergence. In this case our results 
are 77 = 0.977, b t = 34 and S" = 83. For the interpolation 
we have used the values for the velocity and B calculated 
earlier. 

In conclusion, we have found very accurate values for 
the spin velocity, the central charge c and the critical 
exponent 77 for the Heisenberg Hamiltonian with S=3/2 
using the DMRG technique that allowed us to study suf- 
ficiently large chains with negligible error. We found 
c = 1.00 ± 0.01 and r\ = 1.02 ± 0.02 which shows, within 
error bars, that this model belongs to the same universal- 
ity class as the S=l/2 Heisenberg model and the WZW 
model with k = 1. The k = 2S = 3 WZW model with 
c = 9/5 and 77 = 3/5 can be ruled out. These values have 
been obtained by an analysis of the ground-state energy, 
the lowest singlet and triplet excitation energies as well 
as the spin correlation function. We have found an ex- 
cellent agreement between the parameters deduced from 
these different data sets. Apart from a larger coupling 
constant the correlation function of the spin-3/2 chain 
has the same multiplicative logarithmic corrections as for 
spin 1/2. 
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FIG. 4. Scaling of the triplet gap (circles) and the com- 
bined gap A (diamonds). Equations (6) and (8) were used to 
extract the exponent r\ (full lines). 
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